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Abstract 

The chemical equihbration of heavy quarks in a quark-gluon plasma proceeds via annihilation 
or pair creation. For temperatm'es T much below the heavy quark mass M, when kinetically 
equilibrated heavy quarks move very slowly, the annihilation in the colour singlet channel is 
enhanced because the quark and antiquark attract each other which increases their probability 
to meet, whereas the octet contribution is suppressed. This is the so-called Sommerfeld effect. 
It has not been taken into account in previous calculations of the chemical equilibration rate, 
which are therefore incomplete for T < M. We compute the leading-order equilibration 
rate in this regime; there is a large enhancement in the singlet channel, but the rate is 
dominated by the octet channel, and therefore the total effect is small. In the course of the 
computation we demonstrate how operators that represent the annihilation of heavy quarks 
in non-relativistic QCD can be incorporated into the imaginary-time formalism. 
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1. Introduction 



If heavy quarks of mass M in a quark-gluon plasma are initially out of thermal equilibrium, 
they quickly equilibrate kinetically by multiple scatterings with gluons and light quarks [1]- 
[5]. (In this paper we consider an ideal limit in which the plasma lives for a long time.) 
At temperatures T <^ M chemical equilibration is much slower because it requires quark- 
antiquark annihilation or pair creation.^ In fact, given that the probability to find a target 
is Boltzmann suppressed, the chemical equilibration rate Fchem is exponentially small at low 
temperatures, Fchem ~ e~^^'^ [6, 7, 8]. 

Given that the heavy quarks are in kinetic equilibrium they move with non-relativistic 
velocity. When the reacting particles have small relative velocity v, their mutual interac- 
tions can have a large influence on the annihilation or production cross section [9, 10], a 
phenomenon known as the Sommerfeld effect. In perturbation theory this effect would first 
show up in the 1-loop correction to a tree- level cross section ctq: 



For Coulomb-like interactions the Sommerfeld effect manifests itself as a non-vanishing con- 
tribution of 0{as/v). When v becomes as small as as, the 1-loop correction can become 
larger than the tree-level result. The naive loop expansion then breaks down, and the en- 
hanced terms need to be resummed. This effect has to be taken into account for any particle 
reactions close to threshold [11], and has been widely discussed, e.g., in connection with tt 
and gluino and squark pair production in hadronic collisions (see e.g. refs. [12, 13, 14] for 
recent work and references). It may also play an important role in the indirect detection of 
dark matter particles [15]. 

At finite temperature the typical heavy quark velocity is of order v ~ ^/tJM < 1. Thus 
the naive perturbative expansion breaks down for v<as, i.e. T<agM. Then a similar 
rcsummation is needed as in hadronic collisions at zero temperature. The Sommerfeld effect 
in thermal dark matter freeze-out has indeed been discussed in many recent works, such as 
refs. [16, 17, 18], and it may also play a role in certain Icptogcncsis scenarios [19]. However, to 
the best of our knowledge, the Sommerfeld effect has not been taken into account in previous 
calculations of the chemical equilibration rate of heavy quarks. 

In ref. [20], a definition of the heavy quark chemical equilibration rate was given which 
is non-perturbative and thus goes beyond the usual formulation in terms of the Boltzmann 
equation. At leading order it gives the same rate as the Boltzmann equation. To include 
more terms of the perturbative expansion, it would be convenient to use non-relativistic QCD 
(NRQCD [21]) for computing the rate. The pair annihilation of heavy quarks is represented 
in NRQCD by an imaginary part in a coefficient of a 4-fermion operator [22]. Such a complex 

'^Weak interactions are not considered here. 




(1.1) 
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coefficient is related to the analytic structure of a corresponding Green's function. One 
purpose of the present paper is to give a formulation of such operators which can be used in 
the imaginary-time formalism. Subsequently, the NRQCD analysis allows us to disentangle 
the contributions from the colour singlet and octet operators to the heavy quark chemical 
equilibration rate, a necessary first step for discussing the Sommerfeld effect. 

This paper is organized as follows. In sec. 2 we discuss how pair annihilation can be 
incorporated in the imaginary-time formulation of NRQCD, and use this to determine the 
contributions of singlet and octet operators to the chemical equilibration rate. In sec. 3 we 
compute the leading-order chemical equilibration rate taking into account the Sommerfeld 
effect. A brief summary is presented in sec. 4. 



2. Non-relativistic QCD in the imaginary-time formalism 
2.1. General formulation 

NRQCD [21] describes non-relativistic heavy quarks, and gluons and light quarks with mo- 
menta much smaller than M. Therefore the annihilation of a heavy quark-antiquark pair 
cannot be described in terms of the fields of the theory. However, when integrating out the 
scale M, one obtains 4-fermion operators. When this is done in real time, their coefficients 
have an imaginary part which corresponds to the annihilation process [22]. 

The imaginary parts of these coefficients arise from a cut, or discontinuity, of 4-point func- 
tions, viewed as a function of some energy variable u, across the real u-axis. If manipulations 
are carried out in the complex tj-planc, as is necessary e.g. in thermal field theory, then the 
imaginary parts of these coefficients have to be represented in a way which reflects this ana- 
lytic structure. This can be achieved by expressing them in a suitable spectral representation. 

The spectral representation of a 2-point function 11 can be written as 



where the spectral density p is a real and odd function of z. In this form 11 can be evaluated 
both for real lo, which corresponds to real time, and for imaginary lo = iun, where Un 
is a Matsubara frequency, which corresponds to imaginary time. By approaching the real 
w-axis in different ways one obtains different operator orderings (for instance, setting co = 
Rea;-|-iO+ yields a retarded correlator, which for a; S> T is equivalent to the time-ordered one). 
The function n(za;n,k) is purely real. Spectral representations are routinely used in finite 
temperature perturbation theory for Hard Thermal Loop [23, 24] resummed propagators. 




(2.1) 
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The 4-ferinion operators of ref. [22] are 



+HiAoses,)+^^oses,), (2.2) 



Os^So) = ^tr«xx^r«v, 08(^-Si) = v^ar«x-x^^r«v ■ (2.3) 

Here V'jX are 2-component non-relativistic spinors, a are the Pauh matrices, and T" are 
generators of SV{Nc), normalized as Tr [T"T^] = The subscripts 1, 8 refer to singlet and 
octet channels, respectively. The absorptive parts of the coefficients read [22] 



Im/i(^S'o) 



2Nr 



ImhCS,) = 0{at), 



Im/8(^So) = ^ 



^N^^'^s+Oia'), ImfseS,) = ^Tra-^ + Oial), (2.4) 

where Cp = {N^ — l)/2Nc. The corrections of 0{a^) are also known, but not needed here. 

The spectral representation of the most general 4-point function involves three energy 
variables instead of a single one as in eq. (2.1). Fortunately, this complication can be avoided 
for the operators of eq. (2.3) at leading order: only the sum of the energies of the annihilating 

particles appears. This is obvious for the s-channcl annihilation. It is also true for the t and 
u-channel annihilation, because the virtual heavy quark is far off-shell, and effectively leads to 
a point-like interaction of the annihilating pair and the produced two gluons. Consequently, 
the operators can be represented in a form similar to eq. (2.1), 



SS 



M 



II 

/•+°°d£ 

7-00 2^ 



/C2,/C3,/C4 



^,(^,+;c.+/C3+^4)-^ V'.*(/Ci) x.(/C2) xK/Ca) V'„(/C4) 



Prstui^) 



Putsri'^) 



where r, s,t,u contain both spin and colour indices, i enumerates the four cases in eq. (2.2), 
X = (t,x) and /C^ = (fc^,kj). Setting e.g. A;?, M-t-zO"*", ki,k2 0, the absorptive parts 
can be read off: 



Im/W 1 r 
— - — ^ — 
M2 ^ 2 



pS«(2M)-p2,,(-2M)J , 

where a suitable choice of indices is understood. Subsequently, computations can be carried 
out also in the imaginary-time formalism, by including 



(2.6) 



E 



Ix ^ 

J-oo 2^ 



^iiK,+K,+Ks+K,yx 



{KuK2,Ks,K4} 
Prstui^) 



Putsri^) 



(2.7) 

jknl + ikn2 — Z iknl + ikn2 + Zj 

in the Euclidean effective action. Here X = {t,:sc), K = (/;:„, k), = J^^^dr J^, ^{k} = 
^Yl{kn} /k' ^{fcn} denotes a sum over fermionic Matsubara frequencies. 
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(a) (b) 

Figure 1: Feynman diagrams for tire computation of the heavy quark chemical equilibration rate 
within NRQCD. Two ovals connected by a dotted line represent the absorptive part of a 4-fermion 
operator; filled squares represent the heavy quark Hamiltonian; and solid lines represent heavy quark 
propagators. The heavy quarks propagate along the imaginary time direction. 



2.2. Definition of the chemical equihbration rate 

Physically, heavy quark chemical equilibration corresponds to the fact that the energy carried 
by kinetically equilibrated heavy quarks is not conserved because of annihilation or pair 
creation; the chemical equilibration rate is a "transport coefficient" describing the average 
non-conservation. Concretely, the task is to compute the connected correlator [20] 

A(r)^£(//(T,xM0,O))^, 0<T<i, (2.8) 

where H denotes the heavy quark Hamiltonian. In terms of the fields in eq. (2.3), the 
Hamiltonian reads H = M{tlj^ip — x^x) + 0{1/M). We expand to first order in the absorp- 
tive action, eq. (2.7), which is 1/M^-suppressed. After a Fourier transformation, A{ujn) = 
/jj^'^'^dr e*'^"'^A(T), and analytic continuation, p^{uj) = Im A(w„ — >• —i[uj + ^0^]), a coefficient 
denoted by f^j,hem ™ I^*-*] extracted as 

Q,hem = , lim uj^[l + 2Uuj)]p^{u) , (2.9) 

i chom < t'' < UJUY 

where F^j^g^^ ~ e~^^/'^ , cjuv ~ T, and /b denotes the Bose distribution. The chemical equili- 
bration rate then follows from F^j^^^^ = il^^^^/ {2XfM'^), where X{ denotes the heavy quark- 
number susceptibility. 

2.3. Perturbative evaluation of the chemical equilibration rate 

Whereas the formulation of sec. 2.2 is in principle non-perturbative (apart from the fact that 
the matching coefficients in eq. (2.4), reflecting ultraviolet dynamics at the energy scale ~ M, 
need to be computed perturbatively) , we now expand in the gauge coupling as well. The free 
heavy quark propagators have the forms 

(V'r(-f^l)V's(^2))o = — 7— — ^ , (x^(Ki)x,(i^2))o = — T ^ , (2.10) 
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where Ek = M + k'^/2M H , and (K) = 1. Feynman diagrams are illustrated in fig. 1. 

Carrying out Wick contractions and Matsubara sums, performing the analytic continuation, 
and taking the cut, we obtain the spectral function 

Pa(^) = ^ E / 1 [PrlriEk, + + - f^sri-Ek, ' Ek, - Uj)] 

X [U{E,^) + UE,^+E,^+u)][UiEk,+u;)-U{E,J] - (a; ^ -a;)} , (2.11) 
where /f denotes the Fermi distribution. Expanding in a small oj, 

ME,^+uj)-MEk,)'^~fAEk,)[l-fAEk,)] , (2.12) 
and omitting exponentially small terms, the coefficient of eq. (2.9) is readily extracted: 

^^chem = -8M2E / [p'rL{Ek,+EkJ-pi%r{-Ek,-Ek,)]MEk,)MEk,)- (2-13) 

i ^ki,k2 

Subsequently we may count the contractions for the operators in eq. (2.3): 

PrssrPsi'So)] ^ 2N,C, , PrssrPsi' S,)] ^ 6iV,Cp . (2.14) 

Identifying the absorptive coefficients from eq. (2.6) (noting that Ef,^ + Ej^^ 2M because of 
the exponential suppression factors), and inserting their values from eq. (2.4), then leads to 

= 16M2 J^^ MEkJfAEkJ X ^^^^ ( i- + + N,y (2.15) 

singlet, octet, Og 

This is the main information needed in the next section. (For completeness we note that a 
division by 2xfM'^ = SN^M^ ^(^fcj ^^^^^ *° ^chem of eq. (3.15) with S^ = S^ = 1.) 



3. Sommerfeld effect in the chemical equilibration rate 

Consider now the annihilation or pair creation of a heavy quark Q and antiquark Q with 
four-momenta /Ci and /C2 . We define v as the velocity of Q in the QQ rest frame: 

_ I I _ ki - k2 
V — \v\ , V = — — — — . (3.1) 
' ' ' 2M ^ ^ 

One has to resum the multiple exchange of gluons with typical momenta Q = {q^,q), where 

~ Mv^ , q=\q\-^ Mv . (3.2) 
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In heavy-quark kinetic equilibrium eq. (3.2) corresponds to 



r , g ~ VmT . (3.3) 



In particular, q is parametrically larger than the Debye scale which is of order gT, where g 
is the gauge coupling (g = ^/A^^as) ■ 

q^gT. (3.4) 

Therefore the Debye screening and Landau damping of the exchanged gluons by the hot 
plasma can be neglected.^ 

The heavy quarks interact with gluons in the plasma, constantly changing their colour 
charge. This could affect the Sommerfeld effect which depends on the colour charge of the 
pair. The scattering with the plasma is characterized by the thermal width 7, which for 
heavy quarks is of order OgT [26]. On the other hand, the virtuality A = (A^ — Q)^ — = 
(A;° — g°)^ — (k — q)^ — of the heavy quark lines is of the same order as the typical 
momentum transfer squared, i.e. A ~ MT. Schematically, a thermal width would replace 

(A:0 _ ^0)2 ^ (^0 _ ^0 + .^^2 ^ (^0 _ ^0^2 ^ 2iA:% (3.5) 

in the propagator. Since A;'^7 ~ UgMT <^ MT ~ A, the width and correspondingly the 
colour change due to scattering with the heat bath are small compared with virtuality, and 

can be neglected at leading order. 

In ref. [20] it was shown that, ignoring the Sommerfeld effect, the leading order chemical 
equilibration rate can be obtained from a Boltzmann equation which contains the Born cross 
section. The rcsummation of the Sommerfeld-enhanced terms modifies the Born matrix 
elements as [9, 10, 11] 

l-A^resummedP = S \Mtieef , (3.6) 



where S = S{v) is the so-called Sommerfeld factor. When the QQ pair is in a colour singlet 



state the Sommerfeld factor is S = S-^ with 



whereas for the octet S = Sg with 



^This is true not only parametrically but also numerically: we have checked that, above threshold and for 
typical parameter values, eq. (3.7) is in excellent agreement with the ratio of resummed and tree-level singlet 
spectral densities [25], in which the effects of Debye screening and Landau damping are included. 
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a, sqrt(M/T) 



1.0 



Figure 2: The averaged Sommerfeld factors, eq. (3.16), for the singlet and octet contributions. 



At tree level the processes gg o QQ and qq o QQ contribute to the chemical equilibration 
rate. The result of ref. [20] can be written as 



tree 

gg^QQ 



(3.9) 



where the sums are over all spin and colour degrees of freedom. In qq -f-)- QQ the QQ is in 
a colour octet state, whereas the process gg o- QQ has both octet and singlet contributions. 
Denoting by r the ratio of octet to singlet contributions, eq. (2.15) implies that 

_ -4. _ 5 
^ ~ 2 ~ 2 ■ 



(3.10) 



According to eq. (3.6) one has to replace in eq. (3.9) 



A^tree _ 


2 

— )- 


A^trec _ 


A^tree _ 


2 

— )- 


A^tree _ 



2/1 r 

Si + /Ss 



1 + r 



1 + r 



The summed tree-level matrix elements in the non-relativistic limit are 

2 



Eh 
El-^ 



tree 
gg^QQ 

tree 



4(7*C7pAr,(4Cp - N,) , 
4/C7p7V, . 



(3.11) 
(3.12) 

(3.13) 
(3.14) 
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Thus we find 



-M/T 



87rM2 V 27r 

with the thermally averaged Sommerfeld factors 



5. = -1 (^) j^v ^^e-^^VT a G {1, 8} . 



, (3.15) 



(3.16) 



After a rescaling of f, the Sommerfeld factors of cq. (3.16) are seen to be functions of 
g^-^M/T only. A numerical evaluation is shown in fig. 2. Analytically, for T <^ a^M we get 



^1 




(3.17) 



whereas Sg is exponentially small (although decreasing only slowly in fig. 2). For T S> a^M, 
on the other hand, 



5i « 1 + 



s 



1 



g\N, - 2Cp 




(3.18) 



As an example, if we take Og ~ 0.3, M ~ 1.5 GeV, and T 2± 300 MeV, then 2± 3.4, 
Sg 2± 0.8. For N{ = 3, this implies that the square brackets in eq. (3.15) evaluate to 4.28 
rather than the naive 4.17. In other words, the substantial Sommerfeld enhancement of the 
singlet channel is all but compensated for by the fact that most channels, in particular all 
associated with light quarks, are octets, and for octets there is a mild suppression. 



4. Summary 

In a heavy ion collision, the heavy quark chemical equilibration rate parametrizes the rate at 
which heavy quarks and antiquarks, produced in overabundance in an initial hard process, 
annihilate during the thermal stage of the fireball evolution. It can be viewed as a fundamental 
property of thermal QCD, whose systematic understanding may have interesting theoretical 
relations to cosmology, given that similar (co-) annihilation phenomena lie e.g. at the heart 
of computations determining the dark matter relic abundance (in some scenarios). 

On general grounds, the pcrturbativc expansion for the chemical equilibration rate has the 
same functional form as the cross section shown in eq. (1.1) with v ~ y/T/M. In this paper, 
we have resummed the terms of 0(a'^/v"'), describing the Sommerfeld effect, to all orders. 
The result has the form shown in eq. (3.15), with numerical factors plotted in fig. 2. Due to a 
fortuitous cancellation between a strongly enhanced but mildly weighted singlet contribution, 
and a mildly suppressed but strongly weighted octet contribution, the numerical results turn 
out to be largely insensitive to the resummation. 
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The cancellation is peculiar to Nc = 3. For instance, for the fundamental representation 
of SU(2), possibly relevant for dark matter (co-)annihilation at temperatures above the elec- 
troweak scale, the repulsive non-singlet contribution is absent (cf. eq. (3.10)). There is only 
an attractive channel also for oppositely charged particles in U(l), and indeed the Sommer- 
feld effect is likely to play an important role in chemical equilibration in hot QED plasmas 
(see e.g. ref. [27] for a general discussion of the problem). 

Even though the 0{as/v) contribution in eq. (1.1) is insignificant in practice for rdigm, the 
functions 0{aslnv) and 0{as) might well be large. Therefore their determination, as well 
as a fully non-perturbative study of the chemical equilibration rate remain, in our opinion, 
well-motivated challenges. 
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